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This  report  describes  an  effort  carried  out  by  the  National  Weather  Service 
(NWS),  with  funding  support  from  the  Federal  Aviation  Administration  (FAA), 
to  develop  an  improved  technique  for  the  prediction  of  ceiling  and  visi¬ 
bility.  The  objective  of  the  effort  was  to: 

(1)  Develop  equations  for  predicting  the  probability  of  occurrence 
of  specific  categories  of  ceiling  and  visibility  for  time  projections 
of  3,  6,  9,  12  and  15  hours  for  a  specified  list  of  air  terminals. 

(2)  Develop  the  necessary  computer  programs  to  operationally 
implement  the  prediction  equations. 

The  probability  forecasts  were  to  be  based  only  on  surface  observational 
data  from  the  terminal  for  which  the  forecast  was  made  and  were  to  be 
produced  within  one  hour  after  observation  time.  The  prediction  equations 
were  to  be  devised  and  tested  with  the  approach  discussed  in  Report  Number 
FAA-RD- 70-26  (Allen,  1970). 

In  developing  the  required  equations,  the  statistical  technique  of  screening 
regression  was  used.  From  a  large  set  of  possible  predictors,  a  subset 
was  chosen  by  a  screening  algorithm  for  a  group  of  predictand  categories. 

The  chosen  predictors  for  each  given  predictand  became  variables  in  a 
mathematical  expression  referred  to  as  a  "single-station"  equation.  For 
each  terminal,  a  unique  set  of  single-station  equations  was  derived. 

The  original  work  along  these  lines  began  with  the  efforts  of  Enger, 
et  al.  (1962),  at  the  Travelers  Research  Center.  Subsequent  investigations 
by  Miller  (1964)  and  Enger,  et  al.  (1964),  demonstrated  the  feasibility 
of  the  method  in  producing  forecasts  of  approximately  the  same  accuracy 
as  manually-prepared  forecasts  and  the  best  objective  techniques  then 
available. 

The  REEP  (Regression  Estimation  of  Event  Probabilities'  procedure  was  used 
to  develop  the  prediction  equations.  In  this  scheme,  the  predictors  and 
predictands  are  binary,  i.e.,  they  can  take  on  values  of  0  or  1  only;  if 
the  value  of  a  particular  element  is  within  a  given  range,  the  corre¬ 
sponding  predictor  is  set  to  1;  otherwise,  it  is  set  to  0.  The  regression 
equations  are  developed  stepwise,  successively  selecting  the  best  predictors 
from  a  large  number  available.  This  screening  of  predictors  continues  until 
some  predetermined  number  of  predictors  is  selected  or  until  none  of  the 
remaining  predictors  would  improve  the  relationship  adequately. 

Allen  (1969)  continued  exploring  the  basic  approach  anJ  designed  an  opera¬ 
tional  test  of  the  method  under  field  conditions.  At  this  stage  in  the 
development,  equations  were  derived  by  screening  possible  predictors  from 
not  only  the  terminal  under  investigation  but  also  from  10  to  14  terminals 
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in  the  surrounding  area;  the  results  were  called  "network"  equations.  Equations 
for  eight  major  terminals  were  developed,  and  forecasts  were  generated  in 
real  time  at  a  central  computer  site  and  provided  to  forecasters  responsible 
for  the  Issuance  of  official,  manually-prepared  forecasts.  In  order  to 
determine  their  contribution  to  forecast  improvement,  the  objective  fore¬ 
casts  were  timed  to  reach  the  field  after  the  subjective  forecasts  had  been 
prepared  but  sufficiently  in  advance  of  filing  deadlines.  The  forecaster 
was  therefore  able  to  examine  the  objective  forecast  and  modify  his  original 
forecast  dependent  upon  the  later  information.  For  a  period  of  seven  months, 
objective  forecasts  were  prepared  every  six  hours  and  transmitted  to  the 
field  forecasters.  Records  were  kept  in  the  field  of  the  original  objective 
and  subjective  forecasts  and  modified  subjective  forecasts.  Comparative 
verification  of  the  three  types  of  forecasts  indicated  that  the  objective 
forecasts  were  of  sane  value  as  guidance  in  preparing  terminal  forecasts, 
but  their  value  was  small  and  Irregular  with  respect  to  the  for*  casting 
of  low  ceiling  and  visibility  conditions  in  difficult  weather  situations. 

It  was  suggested  that  these  results  were  due  mainly  to  the  use  of  simple 
predictors  in  the  equations  and  that  using  more  sophisticated  predictors 
would  lead  to  improved  forecasts.  The  experiment  did  demonstrate,  however, 
the  practicalities  of  computing  and  distributing  automated  forecasts  from 
a  central  location. 

Because  of  the  relatively  high  data  and  computer  costs  involved  in  develop¬ 
ing  network  equations,  Allen  (1970)  conducted  an  experiment  in  which  equations 
were  developed  by  screening  possible  predictors  derived  from  variables 
observed  at  only  the  terminal  in  question.  These  single-station  equations 
were  developed  and  tested  for  four  terminals.  In  developing  the  equations, 
screening  was  performed  on  339  variables  including  both  simple  and  compound 
(Boolean)  predictors.  The  latter  were  composed  of  two  or  more  simple 
predictors  connected  by  the  logical  operators  "And"  and  "Or".  Boolean 
predictors  were  used  to  model  physically  meaningful  relationships  which 
cannot  be  described  by  simple  predictors  alone. 

Allen  compared  the  single-station  equations  for  the  four  terminals  with 
network  equations  developed  for  the  same  terminals  in  the  earlier  experiment. 
The  evaluation  was  achieved  by  comparing  the  respective  reductions  of 
variance  attained  during  the  derivation  of  equations  for  each  system.  The 
results  showed  the  network  equations  to  be  slightly  superior  but  not  enough 
so  to  justify  the  significantly  higher  costs  involved  in  developing  equations 
for  a  large  number  of  terminals.  Based  on  this  comparison  and  the  findings 
of  the  earlier  work,  it  was  decided  to  expand  the  single-station  equation 
project  and  develop  equations  for  a  large  group  of  U.S.  terminals. 

DEVELOPMENT  OF  SINGLE-STATION  EQUATIONS 

Discussions  between  the  FAA  and  the  NWS  led  to  a  decision  to  develop  single¬ 
station  equations  for  (1)  20  terminals  which  were  receiving  official  NWS 
forecasts  (FT's)  and  which  were  being  treated  by  the  NWS  in  a  parallel  effort 
—under  the  same  Interagency  Agreement— to  develop  prediction  equations 
utilizing  parameters  output  by  certain  numerical  models;  the  results  of  this 
latter  effort  were  called  MOS  (Model  Output  Statistics)  equations  (Bocchieri 
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Table  1.  Terminals  Selected  for  the  Development  of  Single-Station  Equations 


(a)  Terminals  in  the  NWS 
FT  Program 


1.  Albany,  N.  Y. 

2.  Atlanta,  Ga. 

3.  Baltimore,  Md. 

4.  Buffalo,  N.  Y. 

5.  Nashville,  Term. 

6.  Boston,  Mass. 

7.  Birmingham,  Ala, 

8.  Cleveland,  Ohio 

9.  Cincinnati,  Ohio 

10.  Washington,  D.C.  (National) 

11.  New  York.  N.Y.  (Kennedy) 

12.  New  Orleans,  La., 

13.  Chicago,  Ill.  (Midway) 

14.  Pittsburgh,  Pa. 

15.  Raleigh-Durham,  N.C. 

16.  Savannah,  Ga. 

17.  St.  Louis,  Mo. 

18.  Louisville,  Ky. 

19.  Tallahassee,  Fla. 

20.  Knoxville,  Tenn. 


(b)  Terminals  not  in  the  NWS 
FT  Program 


1.  Bedford,  Mass. 

2.  Jackson,  Miss.  (Hawkins  Field) 

3.  Greenville,  S.C. 

4.  Middletown,  Pa. 

5.  Moses  Lake,  Wash. 

6.  Spartanburg,  S.C. 

7.  Idaho  Falls,  Id, 


and  Glahn,  1972),  (Bocchieri  et  al.,  1973).  The  intent  was  to  compare  the 
two  prediction  systems,  as  well  as  a  third  system— a  hybrid  which  would  com¬ 
bine  the  best  aspects  of  the  two;  and  (2)  7  terminals  which  were  not  receiving 
FTs.  The  27  terminals  selected  for  equation  development  are  shown  in  Table 
1. 

The  single-station  equations  were  to  be  developed  such  that  forecasts  of 
ceiling  and  visibility  would  be  available  for  the  terminals  listed  in  Table 
1(b)  for  3,  6,  9,  12,  and  15-hour  projections,  and  for  the  terminals  listed 
in  Table  1(a)  for  4,  7,  10,  13,  and  16-hour  projections.  The  projections 
for  the  latter  group  were  selected  to  allow  field  forecasters  sufficient 
time  to  use  the  objective  forecasts  as  guidance  material  in  the  preparation 
of  FTs  and  still  meet  filing  deadlines. 
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It  was  also  necessary  that  improved  computer  programs  be  developed  to 
derive  single-station  equations.  The  programs  used  in  the  previous 
effort  were  appropriate  to  that  experimental  work.  However,  they  were 
written  for  a  now  obsolescent  computer  and  required  extensive  human 
intervention.  Furthermore,  the  original  programs  were  designed  to  handle 
data  that  required  considerable  preliminary  processing.  Therefore,  a  new 
set  of  programs  was  developed  which  performs  the  following  operations: 

(a)  Process  hourly  surface  observations  input  in  the  standard  formal 
of  the  National  Climatic  Center.  A  subprogram  reads  the  data,  checks  f,i 
missing  and  erroneous  entries,  tests  for  meteorological  and  chronological 
consistency,  and  finally  converts  the  data  into  binary  or  "dummy"  form. 
This  procedure  examines  13  elements  of  each  observation  and,  for  each 
element,  determines  if  its  value  satisfies  the  criteria  for  each  of  a 
number  of  dummy  variables.  162  observation  dummies  are  associated  with 
observed  elements  as  shown  in  Table  2. 


Table  2.  Observational  Elements  From  Which  Dummy  Variables  Were  Defined 


Element 

Unit  of  Measurements 

Number  of 

Observation  Dummies 

Ceiling  Height 

feet 

9 

Prevailing  Visibility 

miles 

10 

Wind  Direction 

16  compass  points 

17 

Wind  Speed 

knots 

9 

Weather 

types 

12 

Dry  Bulb  Temperature 

°F 

13 

Dew  Point  Temperature 

°F 

5 

Sea  Level  Pressure 

mb 

5 

Total  Cloud  Amount 

tenths 

4 

Relative  Humidity 

% 

6 

Lower  Sky  Cover 

classes 

9 

Time  of  Day 

Local  Standard 

31 

Day  of  Year 

Julian  days 

32 

Total 

162 

(b)  Transform  the  observation  dummy  variables  into  417  "event"  dummy 
variables.  This  procedure  combines  the  162  dummy  variables  for  a  given 
observation  with  the  dummy  variables  derived  from  the  three  preceding 
hourly  observations.  In  this  manner,  an  event  is  created  linking  current 
and  previous  weather  observations. 

(c)  Transform  the  417  event  du.jiies  into  the  329  predictors  which 
are  screened  during  the  derivation  ox  equations.  The  329  predictors  are 
listed  in  appendix  A.  The  transformation  is  accomplished  by  treating 
some  individual  event  dummies  as  predictors  and  by  combining  some  event 
dummies  to  form  compound  (derived)  predictors.  For  example,  predictor 
number  273  is  set  to  1  (true)  if  the  visibility  is  not  currently  greater 
than  4  miles  and  the  wind  direction  is  between  SW  and  WNW  (inch)  and  the 
weather  is  rain  or  rain  showers  or  drizzle  or  freezing  rain. 

(d)  Develop  prediction  equations,  automatically  selecting  the  most 
promising  predictors  from  the  329  available.  A  separate  equation  is 
obtained  for  each  of  five  categories  of  ceiling  and  visibility  (the 
predictands),  to  be  obtained  by  Introducing  each  of  the  available  variables 
into  the  regression.  It  then  selects  the  best  one  and  puts  it  into  the 
regression,  computing  the  appropriate  coefficients.  This  process  is 
iterated  until  30  predictors  are  in  the  regression.  Since  the  input  pre¬ 
dictands  are  either  0  or  1,  the  equations  obtained  give  the  probabilities 
of  occurrence  of  the  several  categories. 

The  computer  program  set  accomplishes  steps  (a)  through  (d)  and  produces 
50  regression  equations  (5  categories  X  5  time  projections  X  2  weather 
elements)  for  a  station;  fifteen  minutes  of  central-processor-unit  (CPU) 
time  on  a  CDC  6600  are  required  if  the  input  record  is  ten  years  long. 

In  operational  practice,  each  of  the  equations  is  evaluated  to  produce  a 
forecast. 


The  form  of  each  equation  is: 


F  -  a  +  a  X.  +  a0  X„  +  ., 
y  o  11  2  2 


.  +  a3Q  X3Q 


where  Fy  is  defined  as  the  probability  of  occurrence  of  the  event  y  given 
the  conditions  represented  by  the  30  selected  predictors  Xj;  aQ  through 
an  are  coefficients.  In  practice,  since  the  predictors  are  binary,  the 
value  of  any  Fy  is  calculated  by  adding  the  coefficients  associated  with 
predictors  whose  value  is  1. 


The  predictand  categories  which  were  chosen  for  this  effort  are  shown  in 
Table  3.  For  each  time  projection  and  for  each  element  to  be  forecast 
(ceiling  and  visibility),  five  equations  were  derived— one  for  each  of  the 
predictand  categories  shown  in  Table  3.  For  a  given  element,  the  predictors 
in  all  five  equations  are  identical  but  the  coefficients  are  unique.  Table 
4  illustrates  this  principle  in  showing  a  portion  of  the  4-hour  ceiling 
prediction  equations  for  Albany,  N.  Y. 
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Table  3.  Definition  of  the  Celling  and  Visibility  Predictand  Categories 


Category 

Ceiling 

(Feet) 

Visibility 

(Miles) 

1 

<  100 

<  3/8 

2 

200  -  400 

1/2  -  7/8 

3 

500  -  900 

1-21/2 

4 

1000  -  1900 

3-4 

5 

>  2000 

>  5 

Table  4.  4-Hour  Celling  Prediction  Equation  for  Albany,  N.Y. 


Predictors 

Predictor  Coefficients 
Predictand  Categories 

for 

1 

2 

3 

4 

5 

Constant 

.001 

.039 

.088 

.202 

.669 

1.  Ceiling  at  t  >  2000  ft 

.001 

-.037 

-.087 

-.196 

.320 

2.  Ceiling  at  tQ  <  400  ft 

.016 

.177 

.009 

-.162 

-.040 

3.  Total  cloud  cover  at  t  >  .9 

o 

-.002 

.004 

.018 

.038 

-.058 

4.  Time  of  day  22-03L  and  RH  at 
t  >  905L  and  weather  at  t 

.054 

.042 

.000 

.005 

-.102 

fog  or  ground  fog  o£  hazeu 
or  smoke 

• 

• 

30.  At  t  ,  ceiling  500-900  ft 
and  wind  direction  NNE-ESE 
and  RH  >  80* 

1 

.040 

-.051 

.000 

-.037 

.047 
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The  final  set  of  329  possible  predictors,  which  was  screened  for  each 
terminal,  is  listed  in  Appendix  A.  This  standard  set  of  predictors  was 
used  to  minimize  the  amount  of  data  processing  required  prior  to  each 
screening  run;  the  predictors  are  applicable  in  any  season,  at  any  time 
of  day,  and 'for  any  terminal. 

The  list  is  composed  of  both  simple  and  Boolean  predictors  and  is  the 
result  of  a  major  effort  to  represent  several  types  of  initial  information. 
With  respect  to  simple  predictors,  it  had  been  found  that  the  initial  value 
of  individual  elements  in  the  observation— especially  the  initial  ceiling 
and  visibility-can  be  a  strong  indicator  of  future  conditions,  at  least 
for  short  time  projections.  On  the  other  hand,  combinations  of  simple 
predictors  were  developed  in  order  to  specify  certain  types  of  initial 
conditions  that  might  presage  low  ceilings  or  low  visibilities,  or  any 
other  range  of  the  predictands.  Some  predictors  were  based  on  changes 
in  ceiling,  visibility,  wind,  or  other  variables  during  the  last  three 
hours  before  forecast  time.  A  large  number  of  predictors  were  Included 
to  represent  different  times  of  day  and  seasons  in  combinations  with  the 
initial  ceiling  or  visibility.  All  parts  of  the  country  were  considered 
in  this  development,  insofar  as  local  forecasting  experience  was  known  or 
had  been  published. 

For  the  derivation  of  equations,  a  data  base  consisting  of  10  years  of 
hourly  surface  observations  (about  88,000)  for  each  terminal  was  used 
whenever  possible;  however,  a  minimum  of  5  years  of  data  was  required 
for  a  terminal  to  be  considered  for  equation  development.  The  most  com¬ 
monly  used  time  period  for  the  data  base  was  from  January  1,  1955  through 
December  31,  1964;  the  period  of  record  of  the  data  base  for  each  terminal 
is  shown  in  Appendix  B.  In  all  cases,  the  data  were  obtained  from  the 
National  Climacic  Center  in  Asheville,  N.C.  and  were  supplied  on  magnetic 
tape  in  the  TDF  1440  format. 

The  computer  programs  which  generate  prediction  equations  from  hourly 
climatic  data  are  available  at  the  NWS. 

MEASURES  OF  PREDICTION  EFFECTIVENESS 

It  was  also  required  that  the  NWS  provide  a  set  of  verification  and 
evaluation  measures  for  the  prediction  equations.  Si’ce  the  equations 
are  to  be  used  to  produce  forecast  guidance  for  airways  forecasters,  the 
effectiveness  of  the  guidance  is  the  extent  to  which  it  improves  the 
weather  forecasts  that  reach  the  aviation  user.  A  direct  measure  of  im¬ 
provement  would  have  required  an  extensive,  controlled  experiment  to 
evaluate  the  difference  between  forecasts  made — under  operational  con- 
ditions--with  and  without  the  guidance.  Such  an  experiment  was  not 
feasible  within  the  resources  available;  however,  by  FAA/NWS  agreement, 
it  was  decided  to  evaluate  several  characteristics  of  the  guidance  that, 
to  a  large  extent,  determine  its  effectiveness: 

1.  Accuracy  3.  Ease  of  Interpretation 

2.  Timeliness  4.  Repeatability 
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Accuracy 


Since  the  forecast  guidance  is  in  terms  of  probabilities,  the  appropriate 
measure  is  the  P-score  (Brier,  1950)1.  is  essentially  the  mean  square 
error  of  a  set  of  probability  forecasts.  Because  it  is  strongly  influenced 
by  the  frequency  of  occurrence  of  the  several  categories  of  ceiling  or 
visibility,  the  score  of  the  climatological  probability  forecast  is  used 
as  a  norm  and  the  accuracy  of  the  forecast  is  expressed  in  percent  improve¬ 
ment  over  climatology.  Thus,  accuracy  for  each  set  of  forecasts  is 
represented  by  the  expression 


P  (dim  )  -  P  (fcst] 
P  (dim  ) 


Timeliness 

It  is  necessary  that  the  guidance  be  in  the  hands  of  the  aviation  forecaster 
when  he  can  make  die  best  use  of  it  in  his  forecast.  In  planning  the 
development  of  the  prediction  equations,  timeliness  of  the  ultimate  guidance 
product  was  considered  essential;  the  guidance  should  be  in  the  hands  of  the 
forecaster  an  hour  before  he  must  issue  his  forecast.  Significant  increase 
of  the  lead  time  would,  of  course,  reduce  the  accuracy  of  the  forecast. 
Although  the  system  is  designed  to  produce  the  guidance  on  schedule,  com¬ 
puter  and  communications  failures  may  cause  the  delay  or  loss  of  some 
guidance  issuances.  Therefore,  timeliness  (and  dependability)  will  be 
measured  by  the  percent  of  guidance  issuances  that  reach  the  forecaster 
one  hour  before  he  must  complete  his  forecast. 

Ease  of  Interpretation 

The  forecast  guidance  product  is  a  statement  of  probabilities  of  the 
several  categories  of  ceiling  and  visibility.  As  with  all  new  forecast 
guidance,  the  field  forecasters  will  be  furnished  background  information 
so  that  they  may  properly  interpret  the  guidance.  In  this  connection, 
explicit  provision  will  be  made  for  inquiries  if  the  forecasters  require 
additional  information.  Also,  a  representative  sample  of  forecasters 
will  be  surveyed  to  determine  whether  the  format  of  the  message  containing 
the  guidance  is  satisfactory.  Ease  of  interpretation  will  be  considered 
satisfactory  if  and  when  all  significant  questions  are  resolved  and  80% 
of  the  forecasters  are  satisfied  with  the  format. 

Repeatability 

Identical  inputs  will  produce  identical  forecast  guidance.  (Note  that 
the  time  and  date  of  the  observations  are  part  of  the  input).  Since 
complete  repeatability  is  inherent  in  the  system,  no  measure  of  repeat¬ 
ability  is  necessary. 

*See  Appendix  F  for  an  explanation  of  all  verification  methods  used 
in  this  report. 
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ANALYSIS  AND  VERIFICATION  OF  THE  EQUATIONS 

Following  their  derivation,  the  equations  for  the  20  terminals  shown  in 
Table  1(a),  were  examined  in  three  ways:  (1)  an  analysis  of  the  predictors 
chosen;  (2)  a  verification  of  the  equations  with  Independent  data;  and 

(3)  an  attempt  to  determine  the  optimum  number  of  predictors  to  be  Included 
in  each  equation. 

The  analysis  of  predictors  was  conducted  to  summarise:  (a)  the  order  of 
selection  of  predictors,  with  respect  to  elements  and  projections,  by  in¬ 
dividual  terminals;  (b)  the  frequency  of  selection  of  each  of  the  329 
possible  predictors  for  the  20  terminals  as  a  group;  and  (c)  the  manner 
in  which  predictors  were  selected  for  elements  and  projections  with  respect 
to  each  of  the  329  possible  predictors. 

Analysis  (a)  consisted  of  a  summary,  by  individual  terminals,  of  the  pre¬ 
dictors  in  each  derived  equation  for  a  given  element  and  projection.  The 
summary,  which  is  presented  in  Appendix  C,  lists  the  30  predictors  Included 
in  each  equation;  the  predictors  are  listed  in  the  order  in  which  they  were 
selected  during  the  screening  process.  For  any  given  element  and  projection, 
the  same  predictors  appear  in  all  5  predictand  category  equations  but  for 
a  given  predictor  the  coefficients  will  vary  from  equation  to  equation. 

Due  to  space  limitations,  the  coefficients  have  been  omitted  from  Appendix  C. 

In  analysis  (b),  the  number  of  times  each  of  the  329  predictors  was  selected 
during  the  derivation  of  equations  was  calculated.  The  maximum  number  of 
times  a  predictor  could  have  been  selected  was  200— if  it  were  included  in 
every  element  and  projection  equation  for  all  20  terminals.  In  actuality, 
the  greatest  number  of  times  any  predictor  was  selected  was  188.  The  most 
frequently  chosen  predictor  was  total  cloud  amount  at  t  (TCA  )  covering 
more  than  9/10  of  the  sky.  The  second  most  frequently  chosen  predictor 
was  present  weather  (WEA  )  observed  as  either  rain  or  freezing  rain  or 
sleet— it  was  selected  1§9  times.  The  third,  fourth,  and  fifth  most  fre¬ 
quently  chosen  predictors— each  was  selected  more  than  100  times— were, 
respectively: 

(3)  Sea-level  pressure  at  tQ  less  than  1010  mb 

(4)  Time  of  day  16-01  LST  and  present  ceiling  below  500  ft 
and  either  present  relative  humidity  greater  than  89%  or 
present  weather  observed  as  drizzle. 

(5)  Time  of  day  22-03  LST  and  present  relative  humidity  greater 
than  89%  and  present  weather  observed  as  either  fog  or  ground 
fog  or  haze  or  smoke. 

A  sumnary  of  analysis  (b)  is  pres'nted  in  two  ways:  (1)  the  number  of  times 
each  predictor  was  selected  appears  in  parentheses  in  Appendix  A  following 
each  predictor  number,  and  (2)  the  100  most  frequently  chosen  predictors 
are  shown  in  Appendix  D. 


In  analysis  (c),  each  of  the  329  possible  predictors  was  examined  to 
determine  the  manner  in  which  it  was  selected— if  at  all— with  respect 
to  both  elements,  all  projections,  and  screening  cycle,  summarized  for  all 
terminals.  Only  predictors  which  were  chosen  at  least  40  times  are  shown 
in  Appendix  E.  The  summary  for  each  predictor  is  arranged  to  permit  rapid 
evaluation  of  its  Importance  at  a  particular  terminal  or  to  a  particular 
projection. 

From  an  independent  set  of  data  forecasts  for  4,  10  and  16  hours  were 
generated  with  the  equations  for  the  terminals  in  Table  1(a)  and  verified 
with  respect  to  three  scores;  the  results  were  then  compared  to  persistence 
and  climatology.  The  data  consisted  of  surface  observations  at  each  terminal 
for  0500,  0600,  0700,  and  0800  GMT  for  each  day  during  the  period  October  1, 
1970  to  March  31,  1971.  For  the  projections  used,  the  forecasts  were  valid 
at  1200,  1800,  and  0000  GMT. 

The  results  of  the  verification,  which  are  shown  in  Table  5,  indicate  that, 
in  general,  the  single-station  equations  produced  the  best  forecasts.  The 
exceptions  were  for  the  4-hour  forecast  where  persistence  scored  the  high¬ 
est  for  percent  correct  (all  categories),  and  the  16-hour  forecast,  where  the 
objective  technique  was  tied  with  climatology  for  percent  correct  (all 
categories).  It  should  be  noted,  however,  that  the  Brier  P-score  is  normally 
the  standard  measure  for  probability  forecasts;  by  this  measure  the  objective 
technique  was  the  best  of  the  three  verified.  In  terms  of  percent  improve¬ 
ment  over  climatology  by  the  forecasts  as  measured  by  the  P-score  (see 
previous  section),  the  objective  technique  again  displayed  skill;  results 
are  presented  in  Table  6.  For  this  verification,  climatology  was  defined 
as  the  relative  frequency  of  each  predictand  category  during  the  one  year 
season  specified  above. 

In  order  to  determ  ne  the  effects  of  reducing  the  number  of  predictors  in 
the  equations,  a  second  test  was  conducted  with  independent  data.  The  test 
consisted  of  generating  forecasts  for  the  winters  (October  1  -  March  31) 
of  1967 '1971  with  equations  containing  from  3  to  30  predictors,  in  multiples 
of  3  predictors.  The  independent  data  sample  was  identical  in  its  makeup 
(i.e.  surface  observations  for  0500-0800  GMT)  to  that  used  in  the  previously 
described  verification;  the  projections  used  and  the  verification  scores 
used  were  likewise  identical. 

The  scores  for  the  comparative  verification  of  ceiling  forecasts  are  shown 
in  Table  7  and  indicate  general  improvement  with  an  increase  in  the  number 
of  predictors;  it  can  be  seen,  however,  that  the  amount  of  improvement 
became  small  when  the  number  of  predictors  approached  30.  The  results  for 
visibility  forecasts  were  similar.  The  conclusion  was  reached  that  using 
more  than  30  predictors  could  not  be  justified,  but,  since  some  improvement 
was  attained  up  to  that  point,  the  number  of  predictors  should  not  be  re¬ 
duced  either. 


Independent  Data  Verification  of  Objective  Forecasts,  Persist* 
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indicates  best  score 


Table  6.  Improvement  Over  Climatology  by  Objective 
Forecasts  with  Independent  Data 


Projection 

Element 

Percent  Improvement 
Over  Climatology 

4-Hr 

Ceiling 

39.2 

Visibility 

31.6 

10-Hr 

Ceiling 

29.0 

Visibility 

19.1 

16-Hr 

Ceiling 

19.9 

Visibility 

16.3 

IMPLEMENTATION  FOR  OPERATIONAL  USE 

In  order  to  make  forecasts  with  the  prediction  equations  for  a  given 
terminal,  the  four  most  recent  hourly  airvays  weather  reports  for  that 
terminal  are  required  as  basic  input.  These  must  be  processed  into 
observation  dummies  and  then  into  event  dummies  and  finally  into  dummy 
predictors  in  the  set  of  329  used  in  the  development.  However,  not  all 
of  the  329  need  be  generated,  but  rather  only  that  subset  of  329  actually 
selected  in  the  KEEP  equations  for  that  station.  The  forecasts  are  then 
obtained  by  summing  those  coefficients  associated  with  predictors  having 
a  value  of  1  (true). 

The  work  involved  requires  the  use  of  a  computer  to  make  the  forecasts 
even  for  a  single  station.  Since  the  large  number  of  constants  involved 
and  the  differences  in  the  equations  from  station  to  station  make  the 
preparation  of  an  implementation  program  for  each  station  a  large  task, 
it  was  decided  to  develop  a  computer  program  which  would  generate  another 
program  or  subroutine  to  make  the  forecasts  for  each  station. 

The  Implementation  program  generator  is  a  program  for  the  CDC  6600  that 
accepts,  as  input,  the  definitions  of  the  various  dummies  and  the  REE? 
equations  for  a  given  station  and  generates  as  output  a  .-omputer  program 
that  will  produce  the  forecasts  for  that  station.  The  first  version  of 
the  program  generator  produced  programs  for  use  on  a  time-shared  computer 
system. 

The  computer  program  to  make  forecasts  for  JFK  (Kennedy  International 
Airport,  New  York)  on  a  time-shared  system  was  implemented  in  the  spring 
of  1971  and  was  made  available  to  forecasters  in  New  York.  However,  since 
the  time-shared  system  had  no  other  access  to  weather  data,  the  individual 
observations  had  to  te  entered  manually.  Although  the  program  minimized 
the  work  required  to  enter  the  data,  the  New  York  forecasters  made  little 
use  of  the  system  and  supported  the  other  approach— computation  of  the 
forecasts  at  the  National  Meteorological  Center  with  automated  distribution 
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Table  7.  Comparative  Verification  of  Celling  Forecasts  Generated 

from  Equations  with  Variable  Numbers  of  Predictors,  Winters 
of  1967-71. 


Number  of 
Predictors 

Verification 

Score 

Brier 

P-Score 

Percent 

Correct 

Allen  Utility 

Score 

(a)  4-Hr  Forecasts  Verifying  at 

1200  GMT 

0 

.3792 

77.82 

437 

3 

.2905 

79.24 

625 

6 

.2849 

79.20 

631 

9 

.2834 

79.33 

631 

12 

.2320 

79.39 

634 

15 

.2796 

79.63 

641 

18 

.2783 

79.73 

643 

21 

.2775 

79.79 

646 

24 

.2760 

79.77 

646 

27 

.2754 

79.88 

643 

30 

.2751 

79.77 

649 

(b)  10-Hr  Forecasts  Verifying  at  1800 

GMT 

0 

.3185 

81.71 

454 

3 

.2687 

81.82 

553 

6 

.2656 

81.91 

559 

9 

.2645 

81.91 

557 

12 

.2632 

81.89 

559 

15 

.2623 

81.94 

560 

18 

.2620 

81.97 

561 

21 

.2615 

81.98 

563 

24 

.2607 

81.96 

563 

27 

.2604 

81.97 

564 

30 

.2596 

82.06 

567 

(c)  16-Hr  Forecasts  Verifying  a 

t  0000  GMT 

1 

0 

.2578 

85.72 

441 

3 

.2353 

85.72 

500 

6 

.2335 

85.74 

505 

9 

.2324 

85.74 

508 

12 

.2324 

85.74 

508 

15 

.2310 

85.75 

511 

18 

.2306 

85.73 

511 

21 

.2304 

85.73 

512 

24 

.2302 

85.72 

512 

27 

.2301 

85.74 

513 

30 

.2300 

85.70 

514 
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to  the  field 


A  second  version  of  the  implementation  £"ierator  was  developed  by  modifying 
the  first.  This  verslou  produces  programs  that  can  be  run  on  an  IBM  360 
computer.  A  separate  subprogram  applies  to  each  station  and  a  main  program 
fetches  the  required  observations  from  the  National  Meteorological  Center's 
IBM  360  data  bank  and  calls  on  the  subprograms  to  make  the  forecasts  from 
the  appropriate  set  of  observations.  The  National  Weather  Service  will 
issue  single-station  forecasts  for  the  20  stations  in  Table  1(a)  on  an 
experimental  basis  for  approximately  six  months,  beginning  in  April  1973. 


DEVELOPMENT  OP  EQUATIONS  FOR  ALASKAN  TERMINALS 


In  the  latter  part  of  1971,  the  Alaska  Region  of  the  National  Weather 
Service  became  interested  in  the  use  of  the  single-station  technique 
for  terminals  in  their  area.  The  Alaska  area  has  some  unique  forecast¬ 
ing  problem^  and,  since  many  techniques  u^ed  in  the  contiguous  U.S. 
which  employ  the  output  of  numerical  models  are  not  adaptable  for  use 
in  Alaska,  the  single-station  method  appeared  to  be  particularly  suitable. 


The  FAA  approved  the  development  of  single-station  equations  for  the  23 
terminals  in  Alaska  shown  in  Table  8.  A  major  difference  from  the  earlier 
work  was  that  one  more  category  was  added  for  both  ceiling  and  visibility. 
The  category  modification  did  not  affect  the  lowest  A  categories  but 
resulted  in  the  following  upper  2  categories: 


Category 

5 

6 


Celling 

(feet) 

2000-4900 

>5000 


Visibility 

(miles) 

5-6 

>7 


The  change  was  requested  by  the  Alaska  Region  to  make  the  forecasts  more 
useful  in  the  mountainous  regions  of  Alaska. 


The  equations  are  to  be  implemented  operationally  by  the  Alaska  Region, 
and  the  forecasts  will  be  routinely  available  at  the  forecast  center  in 
Anchorage  for  guidance  and  other  applications. 


Table  8.  Alaskan  Terminals  for  Which  Single-Station 
Equations  were  Developed 


1.  Anchorage 

9.  Cordova 

17.  Betties 

2.  Fairbanks 

10.  Bethel 

18.  Kenai 

3.  Juneau 

11.  Nome 

19.  Sheyma 

4.  King  Salmon 

12.  Kotzebue 

20.  Kodiak 

5.  Annette 

13.  McGrath 

21.  Barter  Island 

6.  Cold  Bay 

14 .  Barrow 

22.  Adak 

7.  Sitka 

15.  Unalakleet 

23.  Summit 

8.  Yakutat 

16.  Northway 
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SUMMARY 


Multiple  linear  regression  equations  were  derived  ror  predicting  the 
probability  of  specified  celling  and  visibility  categories  at  50  terminals. 
The  equations  are  based  upon  weather  observations  at  the  local  terminal 
only  and  were  derived  by  using  the  REEP  screening  technique  on  329  possible 
predictors  consisting  of  simple  and  Boolean  types*  The  data  base  for 
screening  was  generated  from  5  to  10  years  of  hourly  observations  or  each 
terminal.  A  special  computer  program  was  developed  which  generates  the 
data  base  and  derives  the  prediction  equations  in  15  minutes  of  computer 
CPU  time  for  each  terminal. 

The  resulting  equations  for  20  terminals  were  analyzed  to  determine  the 
order,  frequency,  and  manner  in  which  predictors  were  chosen  with  respect 
to  terminal,  meteorological  element,  and  time  projection.  Forecasts  for 
the  same  group  of  20  terminals  were  generated  and  verified  with  one  winter 
season  of  independent  data.  In  general,  forecasts  from  the  single-station 
equations  were  superior  to  both  persistence  and  climatology. 

Equations  for  the  20  terminals  were  also  evaluated  to  determine  the  effects 
of  varying  the  number  of  predictors  in  each  equation.  Comparative  verifi¬ 
cation  of  forecasts  made  from  equations  containing  from  3  to  30  predictors— 
in  multiples  of  3— indicated  that  general  improvement  resulted  by  increasing 
the  number  of  predictors,  but  that  the  amount  of  improvement  became  small 
as  the  predictor  maximum  was  approached. 

A  special  computer  program  was  developed  to  implement— for  operational  use- 
equations  for  about  20  terminals  at  the  National  Meteorological  Center  in 
Suitland,  Maryland.  Equations  for  23  Alaskan  terminals  were  made  available 
to  the  National  Weather  Service  Alaska  Region  for  operational  implementation 
at  the  forecast  center  in  Anchorage. 

CONCLUSIONS 

REEP  prediction  equations,  which  require  only  locally-available  meteoro¬ 
logical  information  to  produce  ceiling  and  visibility  forecasts,  can  be 
dtrived  with  moderate  canputer  and  data  costs.  The  equations  are  adaptable 
to  completely  automated  procedures  and  can  be  used  to  yield  terminal 
forecasts  whenever  four  consu.'utive  hourly  observations  are  available. 

Verification  of  the  prediction  equations  demonstrated  the  soundness  of  the 
single-station  approach.  It  should  be  noted,  however,  that  the  magnitude 
of  improvement  over  persistence  and  climatology  was  not  overwhelming.  This 
result  was  not  unexpected  inasmuch  as  the  single-station  technique  is 
basically  a  form  of  conditional  climatology  and  is  therefore  inherently 
limited  in  skill.  This  research  and  development  technique  provides  an 
Interim  method  for  the  prediction  of  ceiling  and  visibility  until  improve! 
methods,  e.g.,  outputs  of  numerical  weather  models,  are  available.  It  is 
expected  that  equations  combining  both  numerical  and  single-station  methods 
will  produce  forecasts  of  superior  accuracy;  indeed,  this  has  already  been 
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Indicated  In  a  limited  test  by  Bocchierl  et  al.  (1973).  Nevertheless, 
the  single-station  prediction  technique  will  continue  to  have  applicability 
in  those  areas  where  numerical  prediction  products  are  not  available  or  in 
operational  situations  where  reliance  on  numerical  products  will  cause 
intolerable  delays  in  disseminating  a  current  forecast  to  the  user. 
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APPENDIX  A 


BINARY  PREDICTORS  SCREENED  FOR  SINCLE-STATION  PREDICTION  EQUATIONS 


INTRODUCTION 

The  329  predictors  listed  in  this  Appendix  are  two-valued  variables.  Each 
variable  has  the  value  1  if  all  the  conditions  specified  in  its  definition 
are  satisfied,  otherwise  it  hrs  the  value  0.  Prediction  equations  are 
derived  by  regression  screening  applied  to  these  predictors. 

DEFINITIONS 

Each  predictor  has  a  serial  number  which  identifies  it  in  the  computer 
screening  runs.  These  identifiers  were  assigned  by  Allen  (1970)  and  all 
those  predictors  listed  by  Allen  were  used  in  this  wo-!:  except  ncs.  18-25, 
283  and  294  which  Involved  time  lags  of  12  or  more  hours. 

The  weather  elements  from  which  the  predictors  are  formed  are  the  following: 
CIG  Ceiling  in  feet  above  ground 

VIS  Prevailing  visibility  in  miles 

WDR  Direction  of  the  surface  wind 
WSD  Speed  of  the  surface  wind  in  knots 
DBT  Dry  bulb  temperature  in  °F 
DPT  Dew  point  temperature  in  °F 
RLH  Relative  humidity  in  % 

SLP  Sea  level  pressure  in  mb 

SCL  Lower  sky  condition,  amount  of  sky  covered  by  the  lowest 
cloud  layer.  See  Table  Al  below  for  the  code. 

TCA  Total  cloud  amount  in  tenths  of  sky  covered 

WEA  Weather  observed  at  the  given  hour,  in  twelve  groups.  See 
Table  A2  below  for  definition  of  the  groups. 

DOY  Day  of  year,  beginning  with  1  for  January  1.  Each  period 

defined  by  DOY  runs  circularly,  thus  DOY  341  -  80  means  the 
period  December  7  to  March  21. 

TOD  Time  of  day.  This  variable  is  the  local  standard  hour  of  the 
latest  observation  used  in  making  the  forecast  (forecast  time). 


Al 


The  numbers  In  parentheses  immediately  following  the  serial  numbers  in- 
dlcate  the  total  number  of  times  each  predictor  was  selected  during  the 
screening  for  prediction  equations  for  the  20  terminals  listed  in 
Table  1(a). 

The  subscript  on  each  predictor  component  indicates  the  number  of  hours 
lag  between  the  time  the  variable  is  observed  and  the  time  of  the  latest 
data  uiied  in  the  forecast.  For  example,  CIG  ,  refers  to  celling  observed 
3  hours  prior  to  forecast  time. 

The  Boolean  operators  used  in  defining  the  predictors  are  and  "+". 

*  is  the  symbol  for  AND 
+  is  the  symbol  for  OR 

Parentheses  define  the  order  in  which  operations  must  be  performed. 


Table  Al.  Code  For  SCL,  Lower  Sky  Condition 


Code  Value 

Condition 

1 

Clear 

2 

Thin  scattered  or  partly  obscured 

3 

Scattered 

4 

Thin  broken 

5 

Broken 

6 

Thin  overcast 

7 

Overcast 

8 

Sky  obscured 

A2 


or  WEA,  Observed  Weather  Condition 


Weather  Indicators 


None 


R-, 

R, 

R+ 

RW- 

,  RW,  RW+ 

L-, 

L, 

L+, 

ZL-, 

ZL, 

ZL+ 

s-, 

s. 

s+. 

SP-, 

SP, 

SP+: 

sw- 

,  SW,  SW+ 

T, 

T+, 

A, 

TOR 

SG,  SG+,  IC 


ZR-,  ZR,  ZR+,  IP-,  IP,  IP+,  IPW-,  IPW,  IPW+ 


F,  IF 


GF 


BS,  BN,  BD,  BY 


H,  K,  D 


SINGLE-STATION  PREDICTORS 


1.  (2)  CIG.  <  100 

2.  (3)  CIG.  200-400 

3.  (6)  CIG„  500-900 

4.  (4)  CIG,  1000-1900 

5.  (71)  CIG„  >  2000 

6.  (15)  CIG,  <  400 

7.  (27)  CIG,  <  900 

8.  (69)  CIG,  <  4900 

9.  (2)  VIS,  <  3/8 

10.  (3)  VIS,  1/2- 7/8 

11.  (4)  VIS,  1-2  1/2 

12.  (3)  VIS,  3-4 

13.  (42)  VIS,  >  5 

14.  (4)  VIS,  <  7/8 

15.  (8)  VIS,  <  1  3/8 

16.  (28)  VIS,  <  1  1/2 

17.  (54)  VIS,  >  7 

26.  (23)  WDR,  CALM 

27.  (7)  WDR,  NNE-E  *  WSDJ  <  5 

28.  (10)  WDR.  NNE-E  *  WSD.  <  9 

29.  (12)  WDR.  NNL-E  *  WSD.  4-19 

30.  (6)  WER.  NNE-E  *  WSD.  >  20 

31.  (3)  WDR.  NNE-E  *  WSD.  >  30 

32.  (25)  WDR.  NB-ESE  *  WSD.  <  9 

33.  (42)  WDR.  NE-ESE  *  WSD.  10-29 

34.  (6)  WDR,  NE-ESE  *  WSD.  >  20 

35.  (8)  WDR.  E-SSE  *  WSD,  <  5 

36.  (91)  WDR.  E-SSE  *  WSD.  4-19 

37.  (1)  WDR.  E-SSE  *  WSD.  >  20 

38.  (19)  WDR.  SSE-SW  *  WSD.  4-9 

39.  (63)  WDR.  SSE-SW  *  WSD.  6-19 

40.  (3)  WDR,  SSE-SW  *  WSD.  >  20 


41.  (4) 

WDR0 

SW-W  *  WSD0 

4-9 

42.  (6) 

WDRo 

SW-W  *  WSD0 

>  10 

43.  (7) 

WDR0 

W-NW  *  WSD0 

4-19 

44.  (15) 

WDRa 

W-NW  *  WSD„ 

>  10 

45.  (4) 

WDR„ 

NW-N  *  WSD0 

6-19 

46.  (36) 

WDR„ 

NW-N  *  WSD„ 

>  10 

47.  (0) 

WDR, 

NW-N  *  WSD„ 

>  20 

48.  (31) 

DBT0 

<  29 

49.  (71) 

»BT„ 

30-44 

50.  (25) 

DBTa 

45-64 

51.  (10) 

dbt„ 

>  85 

52.  (0) 

DBT„ 

>  90 

53.  (18) 

DPT0 

<  29 

54.  (4) 

DPT„ 

30-39 

55.  (36) 

DPT„ 

>  60 

56.  (22) 

DPTn 

>  70 

57.  (20) 

8LH„ 

<  49 

58.  (42) 

RLH„ 

<  69 

59.  (19) 

RLH0 

70-89 

60.  (5) 

RLH„ 

90-94 

61.  (48) 

RLH„ 

>  90 

62.  (10) 

RLHQ 

>  95 

63.  (32) 

SLP„ 

>  1024.95 

64.  (115)  SLP, 

,  <  1009.95 

65.  (24) 

SCLo 

1 

66.  (31) 

SCL0 

2,  3 

67.  (8) 

SCL„ 

4,  5 

68.  (2) 

SCLq 

6,  7 

69.  (4) 

SCL„ 

8 

70.  (0) 

SCL0 

<  3 

71.  (0) 

SO L, 

<  5 

72.  (14)  TCA0 

<  1 

73.  (188)  TGA, 

,  >  9 

74.  (33)  TCA„ 

<  5 

A5 


75.  (27)  WEA,  None 

76.  (5)  WEA,  R  +  RW  +  L  +  ZL 

77.  (2)  WEA„  L  i  -  ZL 

78.  (30)  WEA,  S  +  SW  +  IC  +  SG  +  SQ 

79.  (2)  WEA,  R  *  P 

80.  (2)  WEA,  (L  +  ZL)  *  F 

81.  (18)  WEA,  R  +  ZR  +  E  +  EW 

82.  (129)  WEA#  R  +  L  +  ZL  +  ZR+E  +  EW 

83.  (9)  WEA,  (R  +  L  +  ZL  +  ZR  +  E  +  EW)  *  F 

84.  (17)  WEA„  (S  +  SW  +  SP  +  IC  +  SG  +  SQ)  *  F 

85.  (13)  WEAC  F  +  GF 

86.  (60)  k'RA.  BS  +  BD  +  K  +  H 

87.  (21)  WEA,  RW  +  SW+T  +  A  +  ZR  +  E 

88.  (39)  (DOT  341-80)  *  (TOD  2200-0500)  *  (CIG,  <  100) 

89.  (32)  (DOY  341-80)  *  (TOD  2200-0500)  *  (CIG,  200-400) 

90.  (27)  (DOY  341-80)  *  (TOD  2200-0500)  *  (CIG,  1000-1900) 

91.  (45)  (DOY  341-80)  *  (TOD  2200-0500)  *  (CIG,  >  2000) 

92.  (58)  (DOY  341-80)  *  (TOD  2200-0500)  *  (VIS,  <  3/8) 

93.  (16)  (DOY  341-80)  *  (TDD  2200-0500)  *  (VIS,  1/2-7/8) 

94.  (16)  (DOY  341-80)  *  (TOD  2200-0500)  *  (VIS,  3-4) 

95.  (41)  (DOY  341-80)  *  (TOD  2200-0500)  *  (VIS,  >  5) 

96.  (19)  (DOY  341-80)  *  (TOD  0600-1100)  *  (CIG,  500-900) 

97.  (11)  (DOY  341-80)  *  (TOD  0600-1100)  *  (CIG,  >  2000) 

98.  (14)  (DOY  341-80)  *  (TOD  0600-1100)  *  (VIS,  1-2  1/2) 

99.  (25)  (DOY  341-80)  *  (TOD  0600-1100)  *  (VIS,  >  5) 

100.  (4)  (DOY  341-80)  *  (TOD  0200-0700)  *  (CIG,  1000-2900) 

101.  (12)  (DOY  341-80)  *  (TOD  0200-0700)  *  (VIS,  3-4) 

102.  (34)  (DOY  341-80)  *  (TOD  1200-1700)  *  (CIG,  <  100) 

103.  (61)  (DOY  341-80)  *  (TOD  1200-1700)  *  (CIG,  200-400) 

104.  (8)  (DOY  341-80)  *  (TOD  1200-1700)  *  (CIG,  1000-1900) 

105.  (42)  (DOY  341-80)  *  (TOD  1200-1700)  *  (CIG,  >  2000) 

106.  (34)  (DOY  341-80)  *  (TOD  1200-1700)  *  (VIS,  <  3/8) 

107.  (7)  (DOY  341-80)  *  (TOD  1200-1700)  *  (VIS,  1/2-7/8) 

108.  (27)  (DOY  341-80)  *  (TOD  1200-1700)  *  (VIS,  3-4) 
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109.  (35)  (DOY  341-80)  *  (TOD  1200-1700)  *  (VIS0  >  5) 

110.  (38)  (DOY  341-80)  *  (TOD  1600-2100)  *  (CIG0  <  100) 

111.  (19)  (DOY  341-80)  *  (TOD  1600-2100)  *  (CIG0  200-400) 

112.  (58)  (DOY  341-80)  *  (TOD  1600-2100)  *  (CIG,  1000-1500) 

113.  (44)  (DOY  341-80)  *  (TOD  1600-2100)  *  (VIS,  <  3/8) 

114.  (13)  (DOY  341-80)  *  (TOD  1600-2100)  *  (VIS,  1/2-7/8) 

115.  (67)  (DOY  341-80)  *  (TOD  1600-2100)  *  (VIS0  3-4) 

116.  (2)  (DOY  81-160)  *  (TOD  2200-0500)  *  (CIG0  <  100) 

117.  (3)  (DOY  81-160)  *  (TOD  2200-0500)  *  (CIG0  1000-1900) 

118.  (44)  (DOY  81-160)  *  (TOD  2200-0500)  *  (CIG0  >  2000) 

119.  (2)  (DOY  81-160)  *  (TOD  2200-0500)  *  (VIS0  <  3/8) 

120.  (2)  (DOY  81-160)  *  (TOD  2200-0500)  *  (VIS0  3-4) 

121.  (42)  (DOY  81-160)  *  (TOD  2200-0500)  *  (VIS,  >  5) 

122.  (50)  (DOY  81-160)  *  (TOD  0600-1500)  *  (CIG,  >  2000) 

123.  (47)  (DOY  81-160)  *  (TOD  0600-1500)  *  (VIS,  >  5) 

124.  (11)  (DOY  81-160)  *  (TOD  1600-2100)  *  (CIG,  <  100) 

125.  (9)  (DOY  81-160)  *  (TOD  1600-2100)  *  (CIG,  200-400) 

126.  (8)  (DOY  81-160)  *  (TOD  1600-2100)  *  (CIG,  1000-1900) 

127.  (7)  (DOY  81-160)  *  (TOD  1600-2100)  *  (CIG,  >  2000) 

128.  (7)  (DO!  81-160)  *  (TOD  1600-2100)  *  (VIS,  <  3/8) 

129.  (10)  (DOY  81-160)  *  (TOD  1600-2100)  *  (VIS,  1/2-7/8) 

130.  (6)  (DOY  81-160)  *  (TOD  1600-2100)  *  (VIS,  3-4) 

131.  (10)  (DOY  81-160)  *  (TOD  1600-2100)  *  (VIS,  >  5) 

132.  (6)  (DOY  161-260)  *  (TOD  2200-0300)  *  (CIG,  <  100) 

133.  (0)  (DOY  161-260)  *  (TOD  2200-0300)  *  (CIG,  1000-1900) 

134.  (44)  (DOY  161-260)  *  (TOD  2200-0300)  *  (CIG,  >  2000) 

135.  (9)  (DOY  161-260)  *  (TOD  2200-0300)  *  (VIS,  <  3/8) 

136.  (4)  (DOY  161-260)  *  (TOD  2200-0300)  *  (VIS,  3-4) 

137.  (53)  (DOY  161-260)  *  (TOD  2200-0300)  *  (VIS,  >  5) 

138.  (56)  (DOY  161-260)  *  (TOD  0400-1100)  *  (CIG,  1000-1900) 

139.  (83)  (DOY  161-260)  *  (TCD  0400-1100)  *  (CIG,  >  2000) 

140.  (47)  (DOY  161-260)  *  (TOD  0400-1100)  *  (VIS,  3-4) 

141.  (61)  (DOY  161-260)  *  (TOD  0400-1100)  *  (VIS,  >  5) 

142.  (36)  (DOY  161-260)  *  (TOD  1200-1900)  *  (CIG,  >  2000) 
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143.  (21)  (DOT  161-260)  *  (TOD  1200-1900)  *  (VIS. 

144.  (12)  (DOT  261-340)  *  (TOD  2200-0500)  *  (CIG, 

145.  (47)  (DOT  261-340)  *  (TOD  2200-0500)  *  (CIG. 

146.  (35)  (DOT  261-340)  *  (TOD  2200-0500)  *  (VIS, 

147.  (23)  (DOT  261-340)  *  (TOD  2200-0500)  *  (VIS. 

148.  (50)  (DOT  261-340)  *  (TOD  0600-1700)  *  (GIG, 

149.  (22)  (DOT  261-340)  *  (TOD  0600-1700)  *  (VIS. 

150.  (22)  (DOT  261-340)  *  (TOD  1800-2100)  *  (CIG. 

151.  (6)  (DOT  261-340)  *  (TOD  1800-2100)  *  (CIG. 

152.  (0)  (DOT  261-340)  *  (TOD  1800-2100)  *  (CIG. 

153.  (23)  (DOT  261-340)  *  (TOD  1800-2100)  *  (CIG. 

154.  (12)  (DOT  261-340)  *  (TOD  1800-2100)  *  (VIS. 

155.  (23)  (DOT  261-340)  *  (TOD  1800-2100)  *  (VIS. 

156.  (17)  (DOT  261-340)  *  (TOO  1800-2100)  *  (VIS. 

157.  (9)  (DOT  261-340)  *  (TOD  1800-2100)  *  (VIS. 

158.  (8)  CIG.  <  100  *  VIS.  <  3/8 

159.  (24)  CIG.  «  400  *  VIS.  <  1  3/8 

160.  (7)  CIG.  <  100  *  USD.  <  3 

161.  (7)  CIG.  <  400  *  USD.  <  5 

162.  (0)  CIG.  1000-1900  *  USD.  >  20 

163.  (1)  CIG.  >  2000  *  USD.  >  20 

164.  (2)  CIG.  <  100  *  UDR.  NNE-E  *  USD.  <  9 

165.  (5)  CIG.  <  400  *  UDR.  NNE-E  *  USD.  6-19 

166.  (1)  CIG.  <  900  *  UDR.  NNE-ESE  *  USD.  >  20 

167.  (8)  CIG.  <  400  *  UDR.  NE-ESE  *  USD.  6-19 

168.  (1)  CIG.  <  400  *  UDR.  SSE-SW  *  USD.  >  6 

169.  (6)  CIG.  <  1900  *  UDR,  E-SSE  *  USD.  >  6 

170.  (74)  CIG.  <  4900  *  UDR,  E-SSE  *  USD.  >  6 

171.  (2)  CIG.  <  1900  *  UDR,  S-SU  *  USD.  >  6 

172.  (26)  CIG.  <  1900  *  UDR.  SU-U  *  USD.  >  6 

173.  (2)  CIG.  >  2000  *  UDR.  SU-U  *  USD.  >  10 

174.  (0)  CIG.  >  2000  *  UDR.  W-NW  *  USD.  >  10 

175.  (87)  CIG.  <  1900  *  UDR.  U-N  *  USD.  >  10 

176.  (31)  CIG.  >  2000  *  UDR.  U-N  *  USD.  >  6 


>5) 

<  100) 

>  2000) 

<  3/8) 

>  5) 

>  20C0) 

>5) 

<  100) 
500-900) 
1000- 1900) 

>  2000) 

<  3/8) 

1-2  1/2) 
3-4) 

>5) 
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177.  (2)  VIS.  <  7/8  *  WSD.  <  5 

178.  (0)  VIS.  >  5  *  WSD.  >  20 

179.  (4)  VIS.  <  3/8  *  WDR.  NNK-E  *  WSD.  <  9 

180.  (7)  VIS.  <  7/8  *  WDR.  NNE-ESE  *  WSD.  6-19 

181.  (10)  VIS.  <  7/8  *  WDR.  E-SSE  *  WSD.  6-19 

182.  (1)  VIS.  <21/2*  WDR.  E-SSE  *  WSD.  6-19 

183.  (0)  VIS.  <13/8*  WDR.  SSE-SW  *  WSD.  6-19 

184.  (2)  VIS.  1-4  *  WDR.  SSE-SW  *  WSD.  6-19 

185.  (3)  VIS.  <21/2*  WDR.  SW-W  *  WSD.  6-19 

186.  (1)  VIS.  >  5  *  WDR.  SW-W  *  WSD.  >  10 

187.  (17)  VIS.  <4  *  WDR.  WSW-NW  *  WSD.  6-19 

188.  (1)  VIS.  >  7  *  WDR.  W-NW  *  WSD.  >  10 

189.  (34)  VIS.  >  5  *  WDR.  W-N  *  WSD.  >  6 

190.  (74)  CIG.  <  400  *  DPT.  >  60 

191.  (54)  CIG.  >  1000  *  DPT.  <  59 

192.  (13)  CIG.  >  2000  *  DPT.  <  39 

193.  (10)  VIS.  <  7/8  *  DPT.  >  60 

194.  (20)  VIS.  >  *  DPT.  <  39 

195.  (5)  CIG.  <  100  *  RIfl.  >  90 

196.  (26)  CIG.  <  400  *  RLH.  >  80 

197.  (48)  CIG.  >  2000  *  RLH.  <  69 

198.  (6)  CIG.  500-900  *  RLH.  >  80 

199.  (2)  CIG.  200-400  *  RLH.  >  80 

200.  (8)  VIS.  <  3/8  *  RLH.  >  95 

201.  (3)  VIS.  <  7/8  *  RLH.  >  90 

202.  (1)  VIS.  1/2-7/8  *  RLH.  >  80 

203.  (4)  VIS.  1-2  1/2  *  RLH.  >  80 

204.  (0)  VIS.  3-4  *  RLH.  <  69 

205.  (33)  VIS.  >  5  *  RLH.  <  69 

206.  (9)  CIG.  <  100  *  WEA.  F  +  GF 

207.  (6)  CIG.  200-400  *  WEA.  F  +  GF 

208.  (8)  CIG.  <  100  *  WEA.  S  SW 

209.  (11)  CIG.  200-400  *  WEA.  R  +  RW  +  L  +  S  +  SW 

210.  (5)  CIG.  500-900  *  WEA.  S  +  SW 
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211.  (6)  CIG,  500-900  *  WEA.  R+RW  +  L  +  S  +  SW+ZR 

212.  (2)  CIGe  1000-1900  *  WEA.  R  +  RW  +  ZR 

213.  (2)  CIG.  1000-1900  *  WEA.  S  +  SW 

*  •  2I4*  (15)  CIG.  >  2000  *  WEA.  R  +  RW  +  8  +  SW 

215.  (32)  CIG.  >  2000  *  WEA.  None 

216.  (8)  VIS.  <  3/8  *  WEA.  R  +  L  +  F 

217.  (3)  VIS.  <  3/8  *  WEA.  F  +  <& 

218.  (7)  VIS.  <  3/8  *  WEA.  S  +  SW 

219.  (14)  VIS.  <  7/8  *  WEA.  L  +  F 

220.  (22)  VIS.  1  3/8  *  W2A.  R  +  L  +  S  +  F 

221.  (12)  VIS,  1/2-7/8  *  WEA.  R  +  L  +  S  +  F 

222.  (33)  VIS.  1-2  1/2  *  WEA.  R  +  L  +  ZR  +  F 

223.  (7)  VIS.  1-2  1/2  +  WEA.  S  +  SW 

224.  (7)  VIS.  1-7  1/2  *  WEA.  R  +  RW  +  L  +  S  +  SW+ZR 

225.  (1)  VIS.  3-4  *  WEA.  R  +  RW  +  L  +  ZR 

226.  (15)  VIS.  <  4  *  WEA.  R  +  RW  +  S  +  SW 

227.  (20)  VIS.  <  4  *  WEA.  L  +  F  +  GF 

228.  (1)  VIS.  >  5  *  WEA.  R  +  RW  +  L  +  S  +  SW  +  ZR 

229.  (3)  GIG.  <  100  *  WDR.  NNE-E  *  RLHo  >  90 

230.  (6)  CIG.  200-400  *  WDR.  NNE-E  *  RLH#  >  80 

231.  (15)  CIG.  500-900  *  WDR.  NNE-ESE  *  RLH.  >  80 

232.  (95)  CIG.  200-900  *  WDR.  NNE-ESE  *  RLH.  >  80 

233.  (38)  CIG.  <  900  *  WDR.  NNE-ESE  *  WSD.  <  9  *  RLH.  >  90 

234.  (38)  CIG.  <  i900  *  WDR.  E-SSE  *  RIH.  >  80 

235.  (3)  CIG.  <  1900  *  WDR.  S-SW  *  RUJ.  «;  e9 

236.  (13)  CIG,  >  2000  *  WDR.  S-SW  *  RLH.  <  69 

237.  (0)  CIG.  <  1900  *  WDR,  SSE-SW  *  RLH.  >  80 

238.  (2)  CIG.  <  1900  *  WDR.  WSW-NW  *  RLH.  <  79 

239.  (4)  CIG.  >  2000  *  WDR,  WSW-NW  *  RLH,  <  69 

240.  (2)  VIS.  <  3/8  *  WDR.  NNE-E  *  RLH.  >  90 

241.  (16)  VIS.  <13/8*  WDR.  NNE-ESE  *  RLH.  >  80 

242.  (19)  VIS.  <13/8*  WDR.  E-SSE  *  RLH,  >  90 

243.  (2)  VIS,  1/2-7/8  *  WDR.  E-SSE  *  RLH.  >  80 

244.  (7)  VIS.  1-2  1/2  *  WDR,  E-SSE  *  RLH.  >  80 


A10 


2^5.  (0)  VIS0  <  2  1/2  *  WDR,  SSE-SW  *  RLH0  >  80 

246.  (14)  VIS.  >  5  *  WDR,  SW-W  *  RLHC  <  79 

247.  (13)  VISD  <  4  *  WDR,  W-N  *  RLH,  >  70 

248.  (21)  VISa  >  5  *  WDRa  W-N  *  RLH,  <  69 

249.  (4)  VIS0  >  7  *  WDR,  W-N  *  RLHa  <  49 

250.  (6)  CIG0  <  400  *  WSD,  <  5  *  RLHo  >  90 

251.  (10)  VIS,  <  7/8  *  WSD,  <  5  *  RLHo  >  90 

252.  (4)  CIG,  >  2000  *  VIS,  >  5  *  (WSD,  >  20  +  RLH,  <  49) 

253.  (1)  CIG,  <  100  *  WDR,  NNE-E  *  WEA,  L  +  F 

254.  (12)  CIG,  <  400  *  WDR,  NNE-E  *  WEA,  r  +  Rw  +  ZR 

255.  (5)  CIG,  <  400  *  WDRo  NNE-ESE  *  WEA,  S  +  SW 

256.  (5)  CIG,  500-900  *  WDR,  NNE-ESE  *WEA,  R  +  L  +  S  +  SW+f 

257.  (5)  CIG,  1000-1900  *  WDR,  NNE-ESE  *WEA,  R  +  L  +  S  +  SW  +  F 

258.  (19)  CIG,  <  400  *  WDR,  NE-SE  *WEA,  R  +  L  +  S  +  F 

259.  (9)  CIG,  500-900  *  WDR,  NE-SE  *  WEA,  R  +  L  +  S  +  F 

260.  (7)  CIG,  <  400  *  WDR,  SSE-SW  *  WEA,  R+L+S+F 

261.  (5)  CIG,  500-900  *  WDR,  SSE-SW  *  WEA,  R+L+S+F 

2b2.  (7)  CIG,  <  1900  *  WDR,  SW-NW  *  WEA,  R  +  RW  +  S  +  SW 

263.  (3)  CIG,  >  2000  *  WDR,  SW-1W  *  WEA,  Nonti 

264.  (3)  VIS,  <  3/8  *  WDR,  NNE-E  *  WEA,  L  +  F 

265.  (10)  VIS,  <  7/8  *  WDRq  NNE-ESE  *WEA,  R  +  L  +  ZR  +  F 

266.  (23)  VIS,  <  7/8  *  WDR,  NNE-ESE  *  WEA,  S  +  SW 

267.  (5)  VIS,  1/2- 7/8  *  WDR,  NNE-ESE  *  WEA,  R  +  L  +  S  +  ZR 

268.  (10)  VIS,  1-2  1/2  *  WDR,  NNE-ESE -*  WEA„  R  +  L  +  S  +  ZR 

269.  (94)  VIS,  <  4  *  WDR,  NE-SE  *  WEA,  R+L  +  S  +  ZR 
2/0.  (12)  VIS,  <2  1/2*  WDR,  SSE-SW  *  WEA,  R+L  +  S+ZR 

271.  (5)  VIS,  >  3  *  WDR,  SSE-SW  *  WEA,  R  +  RW  +  S  +  SW 

272.  (18)  VIS,  >  5  *  WDR,  S-W  *  WEA,  None 

273.  (3)  VIS,  <  4  *  WDR,  SW-WNW  *  WEA,  R  +  RW  +  L  +  ZR 

274.  (0)  VIS,  <  4  *  WDR,  SW-WNW  *  WEA,  *  WEA,  S  +  SW 

275.  (6)  VIS,  >  5  *  WDR,  W-N  *  WEA,  None 

276.  (24)  CIG,  <  400  *  WSD,  <5*  WEA,  L  +  F  +  GF 

277.  (3)  CIG,  <  100  *  WSD,  <  5  *  WEA,  L  +  F  +  (!F 

278.  (2)  VIS,  <  3/8  *  WSD,  <  5  *  WEA,  L  +  F+GF 


All 


.{ 


I 


g: 

?l 


jyf 


Ifc 
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27S.  (14)  VIS,  <  7/8  *  USD,  <  5  *  WEA,  L  +  P  +  GP 

280.  (23)  TOD  2200-0500  *  RIH,  >  80  *  ((WEA_2  R  +  RW)  +  (WEA_1  R  +  RW)) 

281.  (2)  CIG,  <  4900  *  USD,  <  9  *  RIi!0  >  90  *  RLH_3  <  79  *  TCA,  >  9 

282.  (0)  CIG,  1000-4900  *  RLH,  <  69  *  WEA,  R  +  L 

284.  (113  TOD  1600-0100  *  CIG,  <  400  *  (Rlii,  >  90  +  WEA,  L) 

285.  (23)  TOD  1800-0300  *  USD,  <  5  RLH,  >  90  *  TCA,  <  1 

286.  (1)  WDR_3  S-SW  *  WDR,  WSW-NW  *  WEA,  R  +  L 

287.  (0)  WDR,  ESE-SSW  *  WSD,  <  5  *  WEA_2  None  *  WEA,  R  +  L 

288.  (10$  TOD  2200-0300  *  RLH,  <  90  *  WEA,  P  +  GF  +  H  +  K 

289.  (27)  WEA_2  R  +  RW  +  L  *  WEA,  F  +  GF 

290.  (1)  VIS,  <  4  *  WDR,  NE-S  *  WSD,  4-9  *  ((DPT  3  <  29 

*  DPT,  30-39)  +  (DPT_3  30-39  *  DPT,  40-49)  +  (DPT_3  40-59 

*  DPT,  60-69)  +  (DPT_3  60-69  *  DPT,  >  70)) 

291.  (7)  VIS,  <  4  *  WDR,  NE-S  *  WSD,  4-9 

*  ((DBT_3  <  29  *  DBT,  35-44)  +  (DBT_3  21-34  +  DBT,  35-54) 

+  (DBT^3  35-44  *  DBT,  45-64)  +  (DBT_3  45-54  *  DBT,  55-74) 

+  (DBT-3  55-64  *  DBT,  65-84))  *  WEA,  R  +  L 

292.  (0)  ((WDR_3  NNW-NNE  *  WDR,  ENE-ESE)  +  (WDR_3  NNE-ENE  *  WDR,  ESE-SSE) 

+  (WDR_3  ENE-E  *  WDR,  SE-S))  *  WSD,  4-9  *  ((DPT_3  5.  29 

*  DPT,  30-39)  +  (DPT_3  30-39  *  DPT,  40-59)  +  (DPT_3  40-59 

*  DPT,  60-69)  +  (DPT_3  60-69  *  DPT,  >  70)) 

*  ((SLP-3  >  1019.95  *  SLP,  1009.95-1019.95) 

+  (SLP_3  1009.95-1019.95  *  SLP,  999.95-1009.95) 

+  (SLP_3  999.95-1009.95  *  SLP,  <  999.95)) 

293.  (2)  ((WDR_3  SSW-SW  *  WDR,  SE-SSE 

+  (WDR_3  S-SSW  *  WDR,  ESE-SE) 

+  (WDR_3  SSE-S  *  WDR,  E-ESE) 

+  (WDR_3  SE-SSE  *  WDR,  ENE-E) 

+  (WDR_3  ESE-SE  *  WDR,  NE-ENE)) 

*  WSD,  4-19  *  RLH,  >  70 

*  ((SLP_3  >  1024.95  *  SLP,  <  1024.95) 

+  (SLP_3>  1X119.95  *  SLP,  1009.95-1019.95) 

+  (SLP_3  999.95-1009.95  *  SLP,  <  999.95)) 
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295, 

.  (0) 

CIG_1  500-900  *  CIG„  <400 

296. 

■  (0) 

CIG_1  10C0-1900  *  0IGo  <  900 

297. 

(2) 

CTG-X  >  2000  *  CIG0  1000-1900 

298, 

(1) 

CIG_3  500-1900  *  CIG0  <  400 

299, 

(0) 

CIG_3  1000-1900  *  CIG0  <  900 

300. 

(39] 

>  CIG_3  >  2000  *  QG,  <  1900 

301. 

(1) 

VISm1  1-4  *  VIS0  <  7/8 

302. 

(3) 

vis-;l  1-4  *  VIS0  <  1  3/8 

3C3. 

(7) 

VIS_3  >  3  *  VIS0  <  1  3/8 

304. 

(1) 

VIS_3  >11/2*  VIS0  <  7/8 

305. 

(12) 

VIS_3  >  5  *  VIS0  <  2  1/2 

306. 

(0) 

VIS_3  >  5  *  VISC  4  *  WSD_3  6-19 

*  WSD„  <  5  *  RLH_3  <  79  *  RLH0 

>  80 

307. 

(2) 

VIS_3  >  3  *  VIS0  <21/2*  WSD_3  ■ 

4-9  *  WSD0  < 

*  RLH_3  <  79  *  RLHo  >  80 

308. 

(0) 

WDR0  SSE-WSW  *  VSD_3  <  9  *  WSD„  > 

.  io 

*  ((SLP_3  >  1009.95  *  SLP0  <  1009.95 

+  (SLP_3  >  999.95  *  SLP0  <  999. 

95)) 

309. 

(34) 

WSD0  <  9  *  RLH0  >  70  *  WEA„  R  +  RW  +  L  +  ZR 

310. 

(11) 

WSD0  <  19  *  RLH0  >  80  *  WEA„  S  +  1 

SW 

311. 

(17) 

CIG0  <  900  *  VIS„  <  2  i/2 

312. 

(17) 

CIG0  <  900  *  WDR0  SE-SSW  *  RLH0  > 

80 

313. 

(9) 

CIG0  <  900  *  WDR„  SSW-WSW  *  RLH„  : 

>  80 

314. 

(8) 

C1G0  <  1900  *  UDR0  WNV-N  *  RLH0  > 

80 

315. 

(71) 

WDR.  NE-SE  *  SCL0  4-7 

316. 

(12) 

WDR0  SE-SW  *  SCL0  4-7 

317. 

(24) 

WDR,,  NW-NNE  *  SCL0  4-7 

318. 

(4) 

CTG0  <  1900  *  WDR0  W-N  *  WEA0  S  + 

SW 

319. 

(3) 

GIG0  <  1900  *  WDR0  W-N  *  WR40  R  + 

RW  +  L  +  ZR 

320. 

(9) 

VIS.  <  4  *  WDR„  W-N  *  WEA„  S  +  SW 

321. 

(0) 

VIS0  <  4  *  WDR,,  W-N  *  WEA0  R  +  RW 

+  L  +  ZR 

322, 

(3) 

CIG_3  <  400  *  CIGo  >  500 

323. 

(1) 

CIG_3  <  400  *  CIG0  >  1000 

324. 

(2) 

VIS_3  <  7/8  *  V1S0  >  1  1/2 
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325.  (8)  VIS_3  <13/8*  VI S0  >  1  1/2 

326.  (2)  VIS_3  <13/8*  VIS,  >  3 

327.  (0)  WDR_3  SSW-W  *  WER.  ESE-S 

*  ((DPT-3  <  29  *  DPT0  >  30)  +  (DPT_3  30-39 

*  DPT.  >  40)  +  0>PT_3  40-59  *  DPT0  >  60) 

+  (DPT_3  60-69  *  DPT„  >70))  *  TCA„  >  6) 

328.  (0)  WDR_3  NW-NNE  *  WDR„  NE-ESE 

*  ((DPT  -  <  29  *  DPTC  >  30)  +  (DPT  ,  30-39 

•  J  — 

*  DPT„  >  40)  +  (DPT-3  40-59  *  DPT0  >  60) 

+  (DPT-3  60-69  *  DPT0  >70))  *  TCA.  >  6 

329.  (0)  WDR_3  NW-N  *  WDR.  NNE-E 

*  WEA  _  None  *  WEA.  R  +  L  +  ZR  +  S 

-2 

330.  (3)  WDR_3  SSE-SW  *  WDR0  NE-ESE 

*  WEA  „  None  *  WEA-  R  +  L  +  ZR  +  S 

-2 

331.  (3)  ((CIG-3  <  400  *  CIG„  >  500)  +  (VIS_3  <  7/8 

*  VIS.>  1  1/2))  *  ((RLH_3  >  90  *  RLH0  <  89) 

+  (WSD_3  <  5  *  WSD.  >  6)) 

332.  (13)  CIG.  <  1900 

333.  (6)  VIS.  <4 

334.  (23)  (DOY  261-80  *  ((TOD  02-05  *  ((VIS_3  >  5 

*  VIS_1  <  4  *  VIS„  <  1  3/8)  +  (VIS-3  1  1/2-6 

*  VIS_1  1-4  *  VIS„  <  2  1/2)  +  VIS0  <  3/8)) 

+  (TOD  22-01  *  ((VIS_3  >  5  *  VIS^  <  2  1/2 

*  VIS„  <  1  3/8)  +  (VIS_3  1-4  *  VIS.  <  1  3/8) 

+  VIS,  <  3/8)))) 

+  (DOY  81-260  *  ((TOD  02-05  *  ((VIS^  <  2  1/2 

*  VIS,  <  1  3/8)  +  VIS„  <  3/8)) 

+  (TOD  22-01  *  ((VIS_1  <13/8*  VIS0  <  1  3/8) 

+  VIS0  <  3/8))))  +  (TOD  18-21  *  VIS0  <  3/8) 

335.  (3)  (DOY  261-80  *  ((TOD  04-09  *  ((VIS_3  >  5 

*  VIS^  <4  *  VIS„  <  1  3/8)  +  (VIS_3  1  1/2-6  *  VIS^  1-4 

*  VIS.  <  2  1/2)  +  VIS.  <  3/8))  +  (TOD  00-05 

*  ((VIS  .  >  5  *  VIS  .  <  2  1/2  *  VIS.  <  1  3/8) 

•  J  —  -I  — 

+  (VIS_3  1-4  *  VIS.  <  1  3/8)  +  VIS.  <  3/8)))) 

A 14 


+  (DOT  81-260  *  ((TOD  04-09  *  ((VIS  <  2  1/2 

*  VIS,  <  1  3/8)  +  VIS0  <  3/8))  +  (TOD  00-05 

*  ((VIS.!  <13/8*  VISa  <  1  3/8)  +  VIS0  <  3/8)))) 
+  (TOD  20-01  *  VIS0  <  3/8) 

336.  (3)  (DOY  261-80  *  ((TOD  02-05  *  ((CIG  3  >  2000 

*  CIG.!  <  1900  *  CIG.  <  400)  +  (CIG.3  1000-4900 

*  CIG.!  500-1900  *  CIG0  <  900)  +  CIG,  <  100)) 

+  (TOD  22-01  *  ((CIG.3  >  2000  *  CIG.!  <  900 

*  CIG0  <  400)  +  (CIG.3  500-1900  *  CIG0  <  400) 

+  CIG.  <  100)))) 

+  (DOY  81-260  *  ((TOD  02-05  *  ((CIG.!  <  900 

*  CIG.  <  400)  +  CIG.  <  100))  +  (TOD  22-01 

*  ((CIG.!  <  400  *  CIG.  <  400)  +  CIG.  <  100)))) 

+  (TOD  18-21  *  CIG.  <  100) 

337.  (28)  (DOT  261-80  *  ((TOD  04-09  *  ((CIG  >  2000 

*  C16.!  <  1900  *  CIG.  <  400)  +  (CIG.3  1000-4900 

*  CIG_!  500-2900  *  CIG.  <  900)  +  CIG.  <  100)) 

-f  (TOD  00-05  *  ((CIG.3  >  2000  *  CIG  <  900 

*  CIG.  <  400)  +  (CIG.3  500-1900  *  CIG.  <  400) 

+  CIG.  <  100)))) 

+  (DOY  81-260  *  ((TOD  04-09  *  ((CIG  <  900 

*  CIG.  <  400)  +  CIG.  <  100))  +  (TOD  00-05 

*  ((CIG.!  <  400  *  CIG.  <  400)  +  CIG.  <  100)))) 

+  (TOD  20-01  *  CIG.  <  100) 

338.  (14)  ((DOY  261-80  *  ((TCD  02-05  *  ((VIS  3  >  5 

*  VIS.!  <  4  *  VIS.  <  1  3/8)  +  (VIS_3  1  1/2-6 

*  VIS.!  1-4  *  VIS.  <  2  1/2)  +  VIS.  <  3/8)) 

+  (TOD  22-01  *  ((VIS.3  >  5  *  VIS.!  <  2  1/2 

*  VIS.  <  1  3/8)  +  VIS_3  1-4  *  VIS.  <  1  3/8 
+  VIS.  <  3/8))))  +  (DOT  81-260 

*  ((TOD  02-05  *  ((VIS_!  —  2  ^  *  VJS.  <  1  3/8) 

+  VIS.  <  3/8))  +  (TOD  22-01  *  ((VIS  <  1  3/8 

*  VIS.  <  1  3/8)  +  VIS.  <  3/8))))  +  (TOD  18-21 

*  VIS.  <  3/8))  *  ((DOY  261-80 


A15 


*  ((TOD  02-05  *  ((CIG^  >  3000  *  CIG^  <  1900 

*  CIG„  <  400)  +  (CIG  ~  1000-4900  *  'JIG  ,  500-2900 

*  CIG,  900)  +  CIG8  <  100))  +  (TOD  22-01 

*  ((CIG-3  >  3000  *  CIG_1  <  900  *  CIG8  <  400) 

+  (CIG_3  500-2900  *  CIG„  <400)  +  C1G„  <  100)))) 

+  (DOT  81-260  *  ((TOD  02-05  *  ((CIG_1  <  900 

*  CIG8  <  400)  +  CIG„  <  100))  +  (TOD  22-01 

*  ((CIG_1  <  400  *  CIG0  <  400)  +  CIG8  <  100)))) 

+  (TOD  18-21  *  CIG0  <100)) 

339.  (7)  ((DOT  261-80  *  ((TOD  04-09  *  ((VIS_3  >  5 

*  VIS^  <  4  *  VISa  <  1  3/8)  +  (VIS-3  1  1/2-6  *  VIS  1-4 

*  VIS8  <  2  1/2)  +  VIS„  <  3/8))  +  (TOD  00-05 

*  ((VIS_3  >  5  *  VIS_1  <21/2*  VIS8  <  1  3/8) 

+  (VIS_3  1-4  *  VIS8  <  1  3/8)  +  VIS8  <  3/8)))) 

+  (DOT  81-260  *  ((TOD  04-09  *  ((VIS^  <  2  1/2 

*  VIS8  <  1  3/8)  +  VIS8  <  3/8))  +  (TOD  00-05 

*  ((VIS-i  <13/8*  VIS8  <  1  3/8)  +  VIS8  <  3/8)))) 

+  (TOD  20-01  *  VIS„  <  3/8)) 

*  ((DOT  261-80  *  ((TOD  04-09  *  ((CIG  „  >  3000 

-3  — 

*  CIG_1  <  1900  *  CIG„  <  400)  +  (CIG  3  1000-4900 

*  CIG_1  500-2900  *  CIG„  <  900)  +  CIG„  <  100)) 

+  (TOD  00-05  *  ((CIG_3  >  3000  *  CIG_1  <  900 

*  CIGc  <  400)  +  (CIG_3  500-2900  *  v;IG„  <  400) 

+  CIG8  <  100))))  +  (DOY  81-260  *  ((TOD  04-09 

*  ((CIG_1  <  900  *  CIG„  <  400)  +  CIG„  <  100)) 

+  (TOD  00-05  *  ((CIG^  <  400  *  CIG„  <  400) 

+  CIG0  <  100))))  +  (TOD  20-01  *  CIG8  <  100)) 
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APPENDIX  B 


PERIODS  OF  RECORD  OF  DEPENDENT  DATA 

This  appendix  lists  the  period  of  record  of  the  data  base  used  to  derive 
prediction  equations  for  each  terminal.  All  dates  shown  are  inclusive. 

STATIONS  IN  CONTERMINOUS  U.  S. 


Terminal 

1.  Albany,  N.  Y. 

2.  Atlanta,  Ga. 

3.  Baltimore,  Md. 

4.  Buffalo,  N.  Y. 

5.  Nashville,  Tenn. 

6.  Boston,  Mass. 

7.  Birmingham,  Ala. 

8.  Cleveland,  Ohio 

9.  Cincinnati,  Ohio 

10.  Washington,  D.  C. 

11.  New  York  (Kennedy),  N.  Y. 

12.  New  Orleans,  La. 

13.  Chicago  (Midway),  Ill. 

14.  Pittsburgh,  Pa. 

15.  Raleigh- Durham,  N.  C. 

16.  Savannah,  Ga. 

17.  St.  Louis,  Mo. 

18.  Louisville,  Ky. 

19.  Tallahassee,  Fla. 

20.  Knoxville,  Tenn. 

21.  Bedford,  Mass. 

22.  Jackson  (Hawkins  Field),  Miss 

23.  Greenville,  S.  C. 

24.  Middleton,  Pa. 

25.  i-Ieses  Lake,  Wash. 

26.  Spartanburg*  S.  C. 

27.  Idaho  Falls,  Id. 


Period  of  Record 


January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

5955 

- 

December 

1964 

January 

1955 

•- 

December 

1964 

January 

1955 

December 

1964 

January 

1935 

•m 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

•m 

December 

1964 

January 

1955 

• 

December . 

.1964 

January 

1955 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

January 

1955 

- 

December 

1964 

July  1953  -  June  1963 

November  1952  -  October 

1962 

January 

1955 

- 

December 

1964 

January  1955 

- 

December 

1964 

January 

1948 

- 

December 

1951 

January 

1960 

«• 

December 

1961 

January  1948 

- 

December 

1954 

B1 


Appendix  B  (Cant'd.) 


Terminal 

1.  Anchorage 

2.  Fairbanks 

3.  Juneau 

4.  King  Salmon 

5.  Annette 

6.  Cold  Bay 

7.  Sitka 

8.  Yakutat 

9.  Cordova 

10.  Bethel 

11.  Nome 

12.  Kotzebue 

13.  McGrath 

14 .  Barrow 

15.  Unalakleet 

16.  Northway 

17*  Betties 

18.  Kenai 

19.  Shemya 

20.  Kodiak 

21.  Barter  Island 

22.  Adak 

23.  Summit 


STATIONS  IN  ALASKA 

Period  of  Record 

January  1955  -  December  1964 

January  1955  -  December  1964 

January  1955  -  December  1964 

January  1956  -  December  1964 

January  1955  -  December  1964 

July  1955  -  December  1964 

January  1954  -  December  1963 

January  1954  -  July  I960 

January  1956  -  December  1964 

January  1955  -  December  1964 

January  1955  -  December  1964 

January  1955  -  December  1964 

January  1955  -  December  1964 

January  1955  -  December  1964 

January  1952  -  December  1961 

January  1952  -  December  1954 
January  1958  -  December  1964 

January  1955  -  December  1964 

January  1955  -  December  1964 

September  1959  -  December  1964 

January'  1955  -  December  1964 

January  1957  *  December  1964 

January  1955  -  December  1964 

January  1951  -  December  I960 
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APPENDIX  C 


ORDER  OF  SELECTION  OF  PREDICTORS  BY  TERMINAL 


This  appendix  lists  the  predictors  in  their  order  of  selection  in  each  of 
the  10  equations  for  the  20  terminals  identified  in  Table  1(a).  The  pro¬ 
jection  times  are  in  hours,  and  the  predictors  are  identified  by  number  as 
given  in  Appendix  A. 
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ELEMENT 


CEILING 

EQUATIONS 


PROJECTION 
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TERN 
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visibility" 
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ECUATIONS 
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APPENDIX  £ 


MOST  FREQUENTLY  CHOSEN  PREDICTORS 


This  appendix  lists,  in  order,  the  100  most  frequently  chosen  predictors 
for  ail  10  equations  for  nhe  20  terminals  identified  in  Table  1(a).  The 
predictors  are  identified  by  number  as  given  in  Appendix  A. 
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APPENDIX  E 


ORDER  OF  SELECTION  OF  PREDICTORS  BY  PREDICTOR 


This  appendix  gives  the  order  of  selection  for  each  predictor  chosen  a 
minimum  of  40  times  (the  first  44  listed  in  Appendix  D)  for  each  equation 
and  each  of  the  20  terminals  listed  in  Table  1(a).  The  predictors  are 
identified  by  number  as  given  in  Appendix  A.  The  terminals  are  coded 
as  follows: 


SAV  Savannah,  Ga. 

MSY  New  Orleans,  La. 

RAL  Raleigh-* Durham,  N.  C. 

DCA  Washington  (National),  D.  C 

ATL  Atlanta,  Ga. 

BHM  Birmingham,  Ala. 

KNO  Knoxville,  Tenn. 

BNA  Nashville,  Tenn. 

STL  St.  Louis,  Mo. 

BUF  Buffalo,  N.  Y. 

ALB  Albany,  N.  Y. 

BOS  Boston,  Mass. 

MDW  Chicago  (Midway),  Ill. 

CLE  Cleveland,  Ohio 

BAL  Baltimore,  Md. 

TAL  Tallahassee,  Fla. 

CVG  Cincinnati,  Ohio 

LOU  Louisville,  Ky. 

JFK  New  York  (Kennedy) ,  IJ.  Y. 
PIT  Pittsburgh,  Pa. 
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APPENDIX  F 


W 


DEFINITION  OF  VERIFICATION  SCORES 

This  appendix  defines  the  Brier  P-score,  Allen  utility  score,  and  percent 
correct. 

It  is  desirable  that  statements  of  the  probability  of  a  weather  event  be 
reliable;  that  is,  over  a  period  of  time  the  event  should  actually  occur 
with  the  frequency  implied  by  the  probability  forecast.  It  is  also  de¬ 
sirable  that  the  probabilities  be  as  close  to  zero  or  to  100  percent  as 
possible  when  the  event  does  not  occur  or  does  occur,  respectively.  The 
Brier  P-score  (P)  (Brier,  1950)  measures  these  two  characteristics  of 
probability  forecasts  and  is  given  by 


)' 


(FI) 


where  on  each  of  N  occasions  an  event  can  happen  in  only  one  of  r  possible 
classes,  and  ffj,  t±2>  **•  ^ir  rePre8ent  the  forecast  probabilities  that 
the  event  will  occur  in  classes  1,  2,  ...»  r,  respectively.  If  the  r 
classes  are  chosen  to  be  mutually  exclusive  and  exhaustive. 


r 

£ 

j=i 


«  1 


(F2) 


for  each  and  every  i  =  1,  2,  ...  N.  E^.  .  takes  the  value  1  or  0  according, 
respectively,  to  whether  the  event  occurred  In  class  j  or  not.  For  perfect 
forecasting,  the  Brier  P-score  will  have  a  value  of  zero  and,  for  the  worst 
possible  forecasting,  a  value  of  two. 

The  percent  correct  (PC)  of  the  total  number  of  forecasts  is  computed  by 


PC  = 


Sum  of  correct  forecasts 
Total  number  of  forecasts 


X  100, 


(F3) 


The  computation  of  a  utility  score  involves  the  use  of  a  utility  matrix. 

A  utility  matrix  essentially  shows  a  series  of  weighting  factors  which  are 
meant  to  represent  the  usefulness  or  utility  of  forecasts  to  the  user. 

The  utility  matrix  used  in  this  study  and  shown  in  Table  FI  was  devised 
by  R.  A.  Allen  after  consultation  with  forecasters  at  several  aviation 
forecast  centers.  The  matrix  may  not  be  too  different  fi-om  that  of  an 


FI 


actual  utility  matrix  of  an  airline,  and  it  has  been  used  by  Enger  et  al 
(1962)  in  evaluating  ceiling  height  forecasts  at  seven  terminals. 


Inspection  of  Table  FI  shows  that  a  correct  forecast  of  category  one  re¬ 
ceives  a  weight  of  1.0,  while  a  correct  forecast  of  category  five  is  given 
a  weight  of  only  .15.  Also,  a  near  miss  such  as  a  forecast  category  one 
when  category  two  is  observed  is  weighted  by  .7.  The  Allen  utility  score 
(AUS)  is  given  by 


AUS  * 


Z  £  %,  V 


where  W  represents  the  weights  shown  in  Table  FI  and  Z  represents  the 
values  in  the  corresponding  boxes  of  a  5  by  5  fore cast -observed  contin¬ 
gency  table.  This  score  therefore  has  the  following  characteristics: 
a)  more  credit  is  given  for  correct  forecasts  of  the  lower,  more  opera¬ 
tionally  significant,  celling  categories  than  for  correct  forecasts  in 
higher  categories,  and  b)  some  credit  is  given  for  near  misses. 


Table  FI.  The  Allen  Utility  Matrix  Used  to  Judge  the  Usefulness  of 
Ceiling  and  Visibility  Forecasts. 


Observed 

Category 


Forecast  Category 


